Detecting Chiral Edge States in the Hofstadter Optical Lattice 
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We propose a realistic scheme to detect topological edge states in an optical lattice subjected 
to a synthetic magnetic field, based on a generalization of Bragg spectroscopy sensitive to angular 
momentum. We demonstrate that using a well-designed laser probe, the Bragg spectra provide 
an unambiguous signature of the topological edge states that establishes their chiral nature. This 
signature is present for a variety of boundaries, from a hard wall to a smooth harmonic potential 
added on top of the optical lattice. Experimentally, the Bragg signal should be very weak. To 
make it detectable, we introduce a "shelving method" , based on Raman transitions, which transfers 
angular momentum and changes the internal atomic state simultaneously. This scheme allows to 
detect the weak signal from the selected edge states on a dark background, and drastically improves 
the detectivity. It also leads to the possibility to directly visualize the topological edge states, using 
in situ imaging, offering a unique and instructive view on topological insulating phases. 

PACS numbers: 37.10.Jk, 03.75.Hh, 05.30.Fk 



a. Introduction Recently, synthetic magnetic fields 
[1] and spin-orbit couplings [2] have been realized for ul- 
tra cold atoms using suitably arranged lasers that cou- 
ple different internal states [3]. This opens the path to 
the experimental investigation of topological phases, such 
as quantum Hall (QH) states, topological insulators and 
superconductors, in a clean and highly controllable envi- 
ronment [4, 5]. Topological phases currently attract the 
attention of the scientific community for their remarkable 
properties, such as dissipationless transport and quan- 
tized conductivities [6, 7]. In this context, the recent 
experimental realization of a staggered magnetic field in 
a 2D optical lattice, exploiting laser-induced gauge po- 
tentials, constitutes an important step in the field [8] (cf. 
also [9, 10]). In the near future, large uniform magnetic 
flux should be reachable using related proposals [9-12], 
allowing optical-lattice experiments to explore the Hofs- 
tadter model [13]. The latter is the simplest tight-binding 
lattice model exhibiting topological transport properties: 
well-separated energy bands are associated to non-trivial 
topological invariants, the Chern numbers, leading to a 
quantized Hall conductivity when the Fermi energy is 
located in the bulk gaps [14, 15]. These transport prop- 
erties are directly related to the existence of chiral edge 
states: while bulk excitations remain inert, these gap- 
less states carry current along the edge of the system. 
According to the bulk-edge correspondence [16, 17], the 
Chern numbers characterizing the bulk bands determine 
the number of edge excitations and their chirality, which 
protects the edge transport against small perturbations. 

In view of the experimental progress [8-10], an impor- 
tant issue is to identify observables that provide unam- 
biguous signatures of topological phases in a cold-atom 
framework [18-27]. This is a crucial topic from an exper- 
imental point of view: Measuring the Hall conductivity is 
more difficult than for solid-state systems due to the ab- 
sence of particle reservoirs coupled to the system. Several 



* Electronic address: ngoldmanATulb.ac.be 



proposals exist for measuring the bulk topological invari- 
ants, e.g. based on spin-resolved time- of- flight [22] or 
density measurements [18, 23]. On the other hand, a di- 
rect detection of the edge states, associated with a proof 
of their chiral nature [20, 21], would witness the non- 
trivial topological order associated to cold-atom QH insu- 
lators. In this Letter, we propose a realistic and efficient 
scheme to probe the chiral edge states of the Hofstadter 
optical lattice, which we believe could be extended to any 
ultracold-atom setup emulating 2D topological phases. 

b. The model We study a two-dimensional fermionic 
gas confined in a square optical lattice and subjected to a 
uniform synthetic magnetic field B = Biz [H, 12]. The 
Hamiltonian is taken to be 

Hn — — T \^ r^ r -^ i27r<3E>m -t - -^ h c 

771, n 

+ Y^ Konf(r) cl^,n^m,n^ (1) 

with Cm,n the field operator defined at lattice site {x = 
ma^y = no), m,n G Z, J the tunneling amplitude, 
and Konf(^) a cylindrically symmetric confining poten- 
tial. Eq. (1) describes non-interacting fermions on a lat- 
tice in the tight-binding regime, subjected to a vector 
potential A = (0, Bx, 0) corresponding to <l> magnetic 
flux quanta per unit cell [13]. 

It is instructive to first discuss the spectrum E{ky) 
obtained by solving the Hamiltonian (1) on an abstract 
cylindrical geometry, without the potential Konf(^) {^y 
denotes the momentum along the closed direction y) . For 
^ = p/q, where p, g G Z, the spectrum can be parti- 
tioned in terms of bulk states and topological edge states, 
as shown in Fig. 1(a). Bulk states exist in the absence of 
boundary and form q well-separated subbands [13]. Be- 
cause of its finite boundaries along x, this system also 
features topological edge states, which propagate along 
the edges of the cylinder [16]. They are located within 
the bulk energy gaps, with a quasi-linear local dispersion 
relation, E /h ~ Veky, where Ve is the group velocity. In 
the spectrum represented in Fig. 1 (a), the two bulk 
gaps are associated to the quantized Hall conductivities 



CTi:/ = ±1 (in units of the conductivity quantum), as they 
host a single edge-state branch per edge, with opposite 
chirahty sign('Ue). In the fohowing, we set <l> = 1/3 and 
choose a Fermi energy Ep = —1.5 J located within the 
first bulk gap, so that the fermionic gas forms a QH in- 
sulator with central density n = l/3a^. 

The clear partition of the single-particle spectrum into 
bulk states and topological edge states still holds in the 
experimental planar geometry and in the presence of 
the confining potential Konf(^), taken here of the form 
Konf(^) = J {"f^ l^f^edgeV ' We illustrate the bulk-edge par- 
tition in Fig. 1 (b), which shows the discrete energy 
spectrum eQ, for an infinite circular wall (7 -^ 00), with 
^edge = 13a: Similarly to the cylindrical case, we obtain 
bulk states within the energy bands and edge states (illus- 
trated in Fig. 1(c)) within the bulk gaps. We stress that 
the number of edge excitation branches (relative to the 
number of physical edges) and their chirality do not de- 
pend on the particular geometry, as they are dictated by 
topological invariants associated to the bulk (see [16, 17] 
and Supplementary Material). We verified that the above 
picture remains valid for finite confinements 7 > 2: In 
agreement with their topological nature, edge states sur- 
vive for a sufficiently weak confining potential that does 
not radically perturb the band structure. 
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FIG. 1: (Color online) (a) Single-particle energy spectrum 
E — E{ky) for the cylindrical lattice subjected to a uniform 
flux $ = 1/3, where ky is the momentum along the closed 
direction. The bulk states are found in the purple bands and 
the Hall conductivity cth = ±1 is indicated for ^f ~ ±1.5 J. 

(b) Discrete energy spectrum for a circular infinite wall with 
^edge = 13a. The Fermi energy Ef = — 1.5J is represented by 
a red dotted line and a labels the single-particle states, i.e. 
ho\ipa) = ea\ipa), where ho is the single-particle Hamiltonian. 

(c) The amplitude \ipe{rn, n)|^ corresponding to the edge state 
at Ce ~ — 1.5J = £^F in (b). (d) The probe shape fL{m,n) = 
fhir) used to detect the edge state in (c). 

c. Angular Momentum Spectroscopy : Our aim is to 
design an experimental probe yielding a clear signature 
from the topological states. We note that a finite confin- 
ing potential Fconf(^) generally leads to non-topological 
edge states [20, 28, 29]. However, QH edge states have 
a crucial property, their chirality (i.e. the sign of their 



angular velocity), that allows to distinguish them from 
the bulk and non-topological states by using a probe sen- 
sitive to angular momentum. Bragg spectroscopy, a form 
of momentum-sensitive light scattering, can provide such 
a probe [20, 21], as we now explain. In its usual imple- 
mentation [30, 31], Bragg spectroscopy probes the linear 
momentum distribution. Here, we propose (a) to use 
a spatial mode carrying angular momentum in order to 
probe the angular momentum distribution, and (b) to 
shape the probing lasers to maximize the probability to 
excite an edge state. We consider two lasers denoted 
1, 2 in high-order Laguerre- Gauss modes with optical an- 
gular momenta /i^2, corresponding to the electric fields 
Ei^2(r) = ^/h^fi^^i^{r)exp{-ili^20 - iuJi^2t), where the 
radial mode functions are fi{r) ex (r/ro)'^'e~^ I'^^o^ and 
(r, 0) are polar coordinates. We assume that the beams 
are set off-resonance from a neighboring atomic transi- 
tion, so that spontaneous emission can be neglected. This 
leads to a scattering Hamiltonian 

i^Bragg(i) = hnY^ilUe-^'^^' + rJe^'^^'Yicp, (2) 
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(3) 



Here, the index q = l2 — li represents the amount of angu- 
lar momentum transferred by the probe (in units of h), 
fux)L = h{uJi —UJ2) is the energy transfer, Vt is the Rabi fre- 
quency characterizing the strength of the atom-light cou- 
pling, and the probe profile is /l(^) = (r/ro)^e~^ ^^°/A/l, 
with L=\li\ + I/2I (cf. Fig. 1 (d)). The operator 4 cre- 
ates a particle in the eigenstate \il)cx) of the unperturbed 
single-particle Hamiltonian, i.e. /^ol'^a) =ea|'0a)- 

Solving the time-dependent problem ^0 + ^Bragg(^) 
to first order, we write the many-body wave function as 
|^(t)) ^ ^o(t)|0) + E(M)^fc^W^"'""'I^O, where |0) = 
n^<£;pCj,|0) denotes the groundstate [34], and where 

\kl) = c|.q|0) denotes an excited state with a single 
fermionic excitation {k > Ep^l < E-p) with energy 
hwki = i^k — q) > 0. With the initial conditions 
(60(0) = 1, bki{0) = 0), the number of scattered particles 
(hereafter referred to as excitation fraction) is given by 



N{q,ujL) = Y,\bki{t)\^ = T,,t, 



(4) 



kd 



where the scattering rate Fgc is given by the Fermi golden 
rule (cf. Supplementary Material for details) 
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with 7TtS^^\uj) = {sin{ujt) / uj)'^ . When the Fermi energy 
Ep is set within a bulk gap, and for small intensities 
hft <^ J, the excitation fraction N{q^ujL) probes the 
dispersion relation eg = ee(M) associated to the gap- 
less edge states I'^g) that lie within this gap, where M 
is a quantum number analogous to angular momentum 
(see Fig. 2 (a)). For an optimized probe shape (obtained 
for ro ~ redge/\/|^|/2), this can be deduced from the 



behavior of the overlap integrals /^^ defined in Eq. 3. 
They are represented in the Uki — q plane in Fig. 2(b) for 
7 = oo. At low frequencies uj^i <C J/h^ we find a con- 
tinuous alignment of resonance peaks ulf ^ Oeq. This 
refiects the linear dispersion relation e/h ~ OqM in the 
vicinity of the Fermi energy, and provides the angular 
velocity 6e ^ —0.07J/h and the chirality (i.e. sign(^e)) 
characterizing the edge states in the lowest bulk gap. We 
find that this result is in perfect agreement with a di- 
rect evaluation of the angular velocity [35] . We empha- 
size that the edge states velocity highly depends on the 
boundary produced by the confinement: Oq significantly 
decreases as the potential Konf(^) is smoothened {e.g. 
Oe ^ -0.02J/h for 7 = 10, 0^ ^ -O.OIJ/H for 7 = 2). 
The absence of substantial response for g' > in Fig. 2 
(b) clearly proves that our setup is effectively sensitive 
to the edge state chirality. Naturally, the signal obtained 
by setting the Fermi energy in the second bulk gap, or by 
reversing the sign of the magnetic fiux <l> ^ — <l>, would 
probe the opposite chirality. 

We obtain the excitation fraction N{q^uoL) at fi- 
nite times through a direct numerical resolution of the 
Schrodinger equation [36]. A typical result is presented 
in Fig. 3(a), for q = ±4, emphasizing the three dis- 
tinct regimes of light scattering: "edge-edge", "bulk- 
edge" and "bulk-bulk". The "edge-edge" regime corre- 
sponds to transitions solely performed between the edge 
states close to ^p* A sharp resonance peak is visible at 
cj£®^ ^ O^q ~ 0.3 J/h for q = —4, and stems from four 
transitions between edge states, as sketched in Fig. 2(a). 
Then, at higher frequencies, ul ~ J/h^ small peaks wit- 
ness allowed transitions between the lowest bulk band 
and the edge states located above Ep. Finally, for 
ujl ~ 2 J/h, many transitions between the two neigh- 
boring bulk bands lead to a wide and fiat signal. This 
bulk-bulk response is significant for both g' = ±4, as a 
consequence of the large density of excited states in this 
frequency range. We stress that a well focused probe al- 
lows to significantly reduce any signal of the bulk. In the 
following, we consider the quantity N{q, uoL)—N{—q, ujl)^ 
which is zero for a system with time-reversal symmetry 
(cf. Fig. 3(b)). We have repeated the calculations for 
several potential shapes, finding no qualitative change 
(cf. Fig. 3(c)). Although it is advantageous to use a 
steep confining potential, the signal from the edge states 
is robust, even in a harmonic trap (7 = 2). 

The dispersion relation being almost linear close to 
Ey^ the number of allowed transitions l^^^) -^ I'^/e) scales 
with the probe parameter q in the "edge-edge" regime (cf. 
Fig. 2(a)). Thus, one observes an increase of the peaks 
for increasing values of 1^1 (cf. Fig. 3(d)). We stress that 
this progression only occurs in the "edge-edge" regime, 
namely when q is chosen such that hw'^^^ ~ Me (7 is smaller 
than the energy difference between E^ and the closest 
bulk band. In the case illustrated in Fig. 3(d), the 
"edge-edge" regime is delimited by \qe.e\ ^ 7. Beyond 
l^e-el, the resonance peak enters the "bulk-edge" regime: 
The excitation fraction N{q,ujL) broadens, N{—q,uoL) is 
no longer negligible, and the linear dispersion relation is 
no longer probed. We thus conclude that a moderate 



value (here \q\ ~ 4) is preferable to keep a narrow peak, 
well separated from the broader "edge-bulk" signal. 
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FIG. 2: (Color online) (a) Sketch of the single-particle en- 
ergy spectrum e^ and the transitions l-^^) -^ \^|:k) probed by 
N{q,(jjL), for (jJl ~ (^ki < J/h. (b) The amplitude l/^J^, as 
a function of the probe parameter q and excitation frequency 
Uki, for $ = 1/3, Ef = -1.5J, L = 13 and ro = 5.1a. The 
confining potential is infinite (7 = 00) and redge = 13a. 
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FIG. 3: (a) Excitation fraction N{q,(jjL) versus probe 
frequency, for an angular momentum transfer q = ±4. 
(b) N{q,ujL)-N{-q,ujL) for $ = and $ = 1/3. (c) 
N{q^ujL) — N{—q^ujL) for several shapes of the trapping po- 
tential Konf (^) = J (^/Tedge)^ and q = —4. (d) Increasing 
then broadening of N{q,(jjL) —N{—q,(jjL), for increasing \q\. 
In all the figures Q = 0.05 J/h, t = 20h/J, L = 13, ro = 5.1a, 
^edge = 13a, $ = 1/3, ^F = — 1.5J and 7 = 00 (except in (c)). 

d. Imaging the edge states With respect to exper- 
imental detection, the Bragg scheme described above 
presents a major drawback. The associated signatures in 
the spatial or momentum densities are small perturba- 
tions on top of the strong "background" of unperturbed 
atoms. For a circular system with Fermi radius Rp^ one 
can expect about A/edge ^ A/hOe ^ RFA/h\ve\ edge 
states, with A the bulk energy gap and Ve the group ve- 
locity (cf. above). Using the parameters from Fig. 3(a), 
one finds A/edge ~ 20 while the total number of atoms 
in the calculation is A^ ~ 200. Scaling to more realistic 
numbers for an experiment (A" ~ 10^), and noting that 
the number of scattered atoms is at most a fraction of 
Aedge^ we conclude that one should be able to detect a 
few tens of atoms at best on top of the signal coming 
from ~ 10^ unperturbed ones: This is a significant ex- 
perimental challenge with present-day technology. One 
possibility to avoid this difficulty is to use an alternate de- 



tection scheme, where the probe also changes the atomic 
internal state. The probe signal can then be measured 
against a dark background (without unperturbed atoms) , 
which allows powerful imaging methods to be used {e.g. 
large aperture microscopy, as recently demonstrated for 
quantum gases in optical lattices [32, 33]) . 

We present in Fig. 4(a) a possible implementation suit- 
able for two-electron atoms with ultra-narrow optical 
transitions, inspired by the electron shelving method. In 
the following, we use the particular case of ^^^Yb atoms 
for illustration (see also Supplementary Material). The 
ground g and excited e states have nuclear spin / = 1/2, 
leading to ground {^|, g^} and excited {e^, e^} manifolds. 
The Zeeman degeneracies are split by a real and relatively 
strong magnetic field B ~ 100 G. The states g^ and e^ 
are initially populated, as laser coupling between these 
two states is used to generate the artificial gauge field 
leading to Eq. (1) (cf. methods described in [11, 12]). In 
order to probe the (+) = {g^^ e^} system, one introduces 
a weaker additional laser, coupling g^ ^ e^. A crucial 
point is to ensure that topological edge states have the 
same structure in the initial and final states. To this 
end, the initially unpopulated states (— ) = {^^,e^} are 
also coupled by a laser generating the same gauge field 
as for (+). After the probe pulse, atoms in the g man- 
ifold are dispatched (possibly detected) using an auxil- 
iary imaging transition g ^ f. Crucially, atoms in the e 
manifold are not in resonance with the imaging light and 
are therefore unaffected. The e^ atoms are subsequently 
brought down to the g manifold using, e.g.^ adiabatic- 
passage techniques, and they are finally detected, with- 
out stray contributions from unperturbed atoms in other 
internal states. 
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FIG. 4: (a) Detection scheme for a two-electron atom, where 
g is the electronic ground state, e a met ast able excited state 
and / a second excited state with short radiative lifetime, 
suitable for imaging. The coupling lasers G induce an artificial 
gauge field and the probe laser P allows to detect edge states 
(cf. text). The circles indicate relative populations in each 
state after the probe pulse, (b) Excitation fraction versus 
probe frequency, for several potential shapes, (c)-(d) Density 
y9(_)(?7i, n) for q — ±4, 7 = 10 and ujl — O.lAJ/h. The probe 
and system parameters in (b)-(d) are the same as in Fig. 3. 



To analyze the effect of this probe, we consider a sim- 
plified level scheme with two internal states only, denoted 
by the indices (±). We suppose that only the (+) sec- 
tor is initially populated. The spatial profile of the cou- 
pling laser is similar to the one used for Bragg excita- 
tions, but now the Pauli principle does not restrict the 
available final states, since the state (— ) is initially unoc- 
cupied. The calculation proceeds as before, with initial 

|0) = ne.<£;^4,+ |0) and final \kl) = 4_Q,+ |0) states. 
We write the coupling to the probe as 



^Shelving(t) =hnY^ ll^cl^_Cp^^e-'^-' + h.C, (6) 

a/3 



where the operator ff^ .^x creates a particle of the (±) 

sector in the eigenstate IV^a), and where I^o has the same 

definition as in Eq. (3), since h-^ = h- = ho. We also 
suppose that hCl ^ J to neglect higher order excitations. 

The excitation fraction N{q^ujL) — N{—q^ujL) is rep- 
resented in Fig. 4(b), showing a clear resonance peak 
at low frequencies Hul <C J. Interestingly, this result 
shows that the low-energy regime is still governed by 
the chiral edge states located in the bulk gap, although 
transitions are now allowed for all the states below Ep, 
including the bulk states (see Supplementary Material). 
Indeed, the signal N{q^ujL) — N^—q^uj^) remains small 
and fiat in the "edge-bulk" region, while the chiral "edge- 
edge" peak stands even clearer than in the Bragg case 
(since more "edge-edge" transitions are allowed between 
states of same chirality). By setting the probe parameters 
(q^ujL) close to a resonance peak, one can now populate 
edge states into the (— ) sector and directly visualize them 
using state-selective imaging. The corresponding den- 
sity p(_)(7Ti,n) = (^(t)|n(_)(7Ti, n)l^(t)) is illustrated in 
Figs. 4(c)- (d) for q = ±4. The clear difference between 
the two images, obtained with different signs of q but 
otherwise identical setups, is a direct proof of the chiral 
nature of the edge excitations populated by the probe. 
Thus, we have demonstrated that our state-dependent 
probe (6), combined with in situ imaging, provides an 
efficient method to identify topological edge states in 
a cold- atom experiment. We believe that this method 
could be generalized to any cold-atom setup subjected to 
synthetic gauge potentials that emulate 2D topological 
phases, even in the presence of interactions or disorder. 
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Appendix A: The bulk-edge correspondence and the 
circular geometry 

In this work, we consider the Hofstadter optical lattice 
described by the Hamiltonian in Eq. (1). In an optical- 
lattice experiment, the atoms are confined by a cylindri- 
cally symmetric confining potential, which we take of the 
form 



Konf(r) = J {r/redgeT' 



(Al) 



where the parameters J, redge and 7 are defined in the 
main text. We are interested in the detection of chi- 
ral edge states, which propagate along the circular edge 
of this system. In the solid-state framework, these chi- 
ral states are responsible for the quantum Hall effect 
[6, 7, 16]. In this appendix, we recall how these topo- 
logical edge states are indeed related to the concept of 
topological invariants, which guarantee their robustness 
against small external perturbations. 

The fundamental concept which relates the chiral edge 
states illustrated in Fig. 1 (c) to topological invariants is 
the so-called bulk-edge correspondence [16, 17]. This the- 
orem stipulates that when specific topological invariants 
associated to the bulk bands are non-zero, the presence 
of edge states at the boundaries is guaranteed and that 
their chirality (i.e. their orientation of propagation along 
the edge) is fixed. Moreover, the edge states energies 
are located within the bulk gaps, as illustrated in Fig. 1 
(a). Here, the topological invariants are defined on the 
first Brillouin zone, or equivalent ly on an abstract two- 
dimensional torus T^, as a result of periodic boundary 
conditions applied to both spatial directions (more pre- 
cisely, the topological invariants are defined on a princi- 
ple fibre bundle P(U(1),T^), which is based on the two- 
torus T^). These topological indices are the Chern num- 
bers Njy^ which are associated to each bulk band Ej,{k) 
through the Thouless-Kohmoto-Nightingale-Nijs expres- 



sion (TKNN) [14] 

N, = ^f {dk^u,{k)\dk^u,{k)) - (k^ ^ ky)dk, (A2) 

where \ujj(k)) is the single-particle eigenstate of the 
Hamiltonian with energy Ej^ik)^ and k — (kx^ky) G T^ 
is the quasi- momentum. When the Fermi energy E^ is 
exactly located in a bulk gap, the Hall conductivity is 
directly related to the Chern numbers, 



(a) 



Ef 



-1.5J 



(b) 



^F = +1.5J 



(^H 



= E ^^ 



(A3) 



Ejy<EF 



which can be directly derived from the Kubo formalism 
[14, 15]. Here the conductivity is expressed in units of 
the conductivity quantum. 

For the specific case studied in this work, where we 
set the magnetic fiux $ = 1/3, the bulk energy spectrum 
splits into three energy bands, which have the associated 
Chern numbers Ni = —1^ N2 = 2 and Ns = —1. There- 
fore, when the Fermi energy lies in the first [resp. second] 
bulk gap, the Hall conductivity corresponds to cth = — 1 
[resp. (Jh = +1], as illustrated in Fig. 1 (a)-(b). These 
results, which are derived from the toroidal geometry, are 
summarized in Table I. 

The bulk-edge correspondence dictates the following 
result: if we solve the same model (1) on an open ge- 
ometry (i.e. a system with boundaries), gapless edge 
states will appear in the bulk gaps. Moreover the num- 
ber of edge-state branches is given by the modulus of the 
Hall conductivity in (A3), and their chirality by sign(cri:/). 
This can be easily visualized by solving the model (1) on 
an abstract cylinder (with Konf = 0). The corresponding 
energy spectrum E{ky) is illustrated in Fig. 1 (a) and is 
described in the main text. Since the cylinder has two 
physical edges, we note that the first bulk gap in Fig. 1 
(a) hosts two edge-state branches with opposite orienta- 
tions (i.e. one for each physical edge) [16]. The presence 
of a single edge excitation (per physical edge) is in agree- 
ment with the fact that \crH\ = |A^i| = 1 when the Fermi 
energy is located in this bulk gap: this is precisely the 
bulk-edge correspondence, which is summarized in Table 

n. 

We stress that the number of edge-state branches (per 
physical edge) is independent of the boundary geometry, 
as it is given by a sum of topological indices (A3) asso- 
ciated to the bulk bands. Therefore, when considering 
the realistic circular geometry produced by the confin- 
ing potential Konf(^) in Eq. (Al), one obtains the same 
edge-state structure propagating along the circular edge 
^ = ^edge as the one obtained from the abstract cylinder 
discussed above: the number of edge- state branches and 
the chirality deduced from them are identical, as these 
properties do not depend on the chosen geometry. How- 
ever, let us comment on the fact that the boundaries 
do affect the dispersion relations, and thus the angular 
velocity, of the edge states (cf. main text). Here, the 
bulk-edge correspondence indicates that the lowest bulk 
gap in Fig. 1(b), which corresponds to the circular ge- 
ometry, hosts a single edge-state branch associated to a 
negative angular velocity, since the corresponding Hall 




0.1 0.2 0.3 

ujki{J/h) 




0-114'/ 



0.1 0.2 0.3 

ujki{J/h) 



FIG. 5: (Color online) The amplitude |/^J^, as a function 
of the probe parameter q and excitation frequency ujkh for 
$ = 1/3, L = 13 and ro = 5.1a. The confining potential 
is infinite (7 = cxo) and redge = 13a. (a) ^f = —1.5 J and 
(b) ^F = +1.5 J. The angular velocity of the edge states 
present in the first [resp. second] bulk gap is Oe ~ ^W lo^ ~ 
-OmJ/h [resp. Oe ^ +0.07 J /h]. Thus, the two bulk gaps 
are associated with opposite chiralities, in agreement with 
the bulk-edge correspondence in Table I. 



conductivity an = — 1 is solely governed by the topologi- 
cal expression (A3). Furthermore, the edge-state branch 
present in the second bulk gap corresponds to the op- 
posite chirality, since cth = +1 when the Fermi energy 
is in the highest gap. These results have been verified 
by directly computing the angular velocity of the edge 
states. 



Oe = {i/h)^iljl{m,n)[hoJ]ilJe{m,n), 



Ji 



= -— 2^{^{'^ + I7 ^) — 0{m^ n))ilj*{m^ n)il)e{m + 1, n) 

+ (^(m — 1, n) — 6{m, n))il)l{m^ n)il)e{m — 1, n) 

{6{m^ n + 1) — 6{m^ n))ip*{m, n)ipe{Tn, n + 1) 
(^(tti, n — 1) — 0{m, n))ipl{m, n)ipe{Tn, n — 1), 



i27r^ri 



i27r^m 



where 6{m,n) = atan2((n — A^/2), (tti— A^/2)) and tti, n = 
1, . . . , A^, and also indirectly through the Bragg signals, 
as illustrated in this Appendix (cf. Fig. 5). 



Appendix B: Angular Momentum Spectroscopy 

The interaction between the Bragg lasers and the 
atoms is described by the Hamiltonian 

J^Bragg(i) = nn{Wge-'"^' + W-ge'"-'), (Bl) 



a/3 
1 



Ia0 = T^ UxrM)M^)fL{r)e' 



(B3) 



Here, the operator c]^ creates a particle in the eigenstate 
I'^a) of the unperturbed single-particle Hamiltonian, i.e. 

^olV^a) = ^alV^a) (cf. main text). Solving the time- 
dependent problem Ho-\-HBragg{t) to first order, we write 



the many-body wave function as 

mt)) = bo{tm + ^b,{t)e-^^^'/^\^,), 

^6o(t)|0) + ^6fc,(t)e-*--^|fe/), (B4) 

where |0) = nzy<£; ^t\^) denotes the groundstate at zero 
temperature, and 

|A:0 = |l...l 1...1 I 0...0 1 0...0), (B5) 

where k > Ep^l < Ep and uj^i = {ck — ^i)/h > 0. Here 
we have restricted the fuh Hilbert space to the subspace 
spanned by the ground state and the excited states that 
are coupled to it to first order in the perturbation (Bl). 
Setting the initial condition (6o(0) = l,bki{0) = 0), 
one finds 

(B6) 
where 5±K) = sin(A±t)/A±, A± = Ki±WL)/2. The 
number of scattered atoms, or excitation fraction, is then 
given by 



N{q,UJL) = Y,\bki{t)\\ 
k,l 

k,l 

(B7) 

In the long-time limit, and neglecting the anti-resonnant 
term (ex 6*^^^^), this yields the standard Fermi golden 
rule 

N{q,0JL) = 27rnH Y. \IlfS'^'H^ki-uJL), (B8) 

k>EF,l<EF 



where S'^^^uj) = {l/7rt){sm{ujt)/ujy 



-> 5{uj). The 



expression (B8) emphasizes the explicit relation between 
the excitation fraction N{q,ujL) and the rates |/j 
presented in Fig. 2 (b). 



q 1 2 
ki\ 



At finite times, it is preferable to evaluate the ex- 
citation fraction through a numerical evaluation of the 
Schrodinger equation 



ih — ; — 

dt 



hn Y^ Wkl-nm{t)bnm{t)e' 



(Ekl-Er,m)t/h 



(B9) 



iujLt 



where Wki-nm{t) = {kl\W{t)\nm) and W{t) = WqC''^^ 
W-qC^^^''^. The many-body wavefunction is still re- 
stricted to the first-order subspace but off-resonant terms 
and deviations from the long-time limit are included. For 
the reasonable finite times and small Rabi frequencies 
Vt <C J /h used in our calculations (cf. Figs. (3)- (4)), we 



find that the excitation fraction obtained from a numer- 
ical evaluation of Eq. (B9) is in perfect agreement with 
Eq. (B7). Note that we use excitation times of several 
h/ J ^ typically, which seem experimentally realistic. This 
is long enough to resolve the edge-edge resonance but 
still too short to neglect the broadening due to the finite 
pulse time (cf. Figs. (3)-(4)(b)). 



Appendix C: The Shelving Method 

The scattering Hamiltonian considered in the "shelving 
method" has the form 



sh^ — iojLt 



ffshelving(t) = hn{Wf'e 



(W^f)^e»'^^*), (CI) 
(C2) 



where the operator c^^rj^\ creates a particle of the (±) 
sector in the eigenstate IV^a), and where /^^ has the same 



definition as in Eq. (B3), since /i(_) = h 



(+) 



Hq. 



In this scheme, we suppose that only the (+) sector 
is initially populated, such that the initial and excited 
states have the following forms 



|0) = |l...l I 0...0)(+)|0...0) 



'(-)' 



(C3) 



\kl) = \l...l 1...1 I 0...0)(^,|0...0 1 0...0)(_,, 



Bf 



(C4) 



where we suppose that Q <^ J/h to neglect higher order 
excitations. We note that k is no longer restricted by the 
Pauli principle, such that uj^i = (e/e — ei)/h may now take 
negative values. We follow the same treatment as for the 
Bragg scheme and we obtain the excitation fraction as 

N{q,,0L) = 2nQh Y^ Y\^'ki\'^^'H^ki-ooL). (C5) 

1<Ef k 

which differs from Eq. (B8) by the fact that the final 
states k are now unrestricted. However, we stress that 
the sum over the initial states, X^k^; ^^ ^^- (^^)' ^^ ^^^^^ 
restricted by the Pauli principle: for ujL<^J/h and when 
£^F = — 1-5J, this allows to probe the edge states that are 
located in the first bulk gap only. This important fact 
leads to the asymmetry highlighted in Figs. 4(c)-(d), 
which demonstrates the specific chirality of these edge 
states (cf. main text). 

We finally stress that the condition 



h 



{-)'- 



--h 



{+)'- 



--hn 



is necessary in order to probe the edge-state structure. 
Indeed, if we consider a simpler scheme in which the 
(— ) sector is no longer subjected to a synthetic gauge 
potential, we find that N{q^uoL) — N{—q^uoL) ^ 0. This 
observation shows that our scheme requires that the 
edge states of the (— ) sector should have the same 



Chern numbers associated to the two lowest bulk bands: 
Hall conductivity in the two bulk gaps (from Kubo formula): 
Chirality of the edge states inside the bulk gaps (circular geometry): 



Ni = -1 

(th{Ey G 1st gap) = — 1 

sign(crjj(lst gap)) = (-) 



A^2 = 2 
gh^Ey G 2nd gap) = +1 
sign(crjf (2nd gap))^(+) 



TABLE I: The Hofstadter model with $ = 1/3: Chern numbers, Hall conductivity and edge-state configurations. 



Two-dimensional Geometries: 


torus (abstract) 


cylinder (abstract) 


circular (realistic) 


Number of ID boundaries: 

Chern number of the bulk band Eiy{k): 

Hall conductivity for ^f G ^th bulk gap (from Kubo formula): 





2 


1 


Number of edge- state branches in the rth bulk gap: 
Chirality of the edge states located in the rth bulk gap: 




2la^"^l 
sign(cr^^) 


\^H \ 

sign(cr}jO 



TABLE H: The bulk-edge correspondence: Topological Chern numbers (defined in the bulk) and the properties of gapless 
edge excitations for the three geometries described in the text. Note that for the cylinder, the edge states at the two different 
physical edges have opposite velocities [16, 17]. 



chirality than the initially populated edge-states of the 
(+) sector, i.e. both systems should be subjected to the 
same magnetic flux (cf. Appendix D). 



Appendix D: Detection scheme using state-changing 
transitions: example for ^^^Yb atoms 

We give here a more detailed account of the detection 
scheme using ^^^Yb atoms (see Fig. 4 (a)). The ground 
g and metastable excited e states have zero electronic an- 
gular momentum but nuclear spin / = 1/2. We denote 
the Zeeman manifolds {^|, qaJ] and {e^, e^} in the ground 
^S'o and ^Pq excited states, respectively. The states g^ 
and e^ are initially populated, as laser coupling between 
these two states is used to generate the artificial gauge 
field [11, 12] leading to Eq. (1). A crucial point is to en- 
sure that topological edge states have the same structure 
in the initial and final states (cf. Appendix C). To this 
end, the initially unpopulated states (— ) = {5'i,e|} are 
also coupled by a laser generating the same gauge field as 
for (+). The degeneracies are split by a relatively strong 
magnetic field, l^Ea^ = —garriiB, where a = e^g denotes 
the ground or excited manifold, z =t / |, m^ = ±1/2 
the nuclear spin quantum number, gg/h ~ —750 Hz/G, 
and ge/h ^ -1250 Hz/G. A bias field P - 100 G thus 
leads to Zeeman shifts ~ ±25 kHz on the tt transitions 
and ^ ±100 kHz on the a^ transitions: these shifts are 
very large compared to the Rabi frequencies of both the 



gauge-field and probe lasers, which can thus be treated 
independently. 

In order to probe the (+) = {^^, e^} system, one intro- 
duces a weaker additional laser, coupling g^ -^ e^. Due 
to the gauge coupling, a population will build up in the g^ 
state as well (roughly equal to that in the e^ state since 
Q <^ J). Those atoms will be missing in the final de- 
tection step. After probing, the lattice sites are isolated 
by rapidly raising the lattice height and switching off the 
artificial gauge field. Atoms in the g manifold are dis- 
patched (possibly detected) using an auxiliary imaging 
transition g ^ f. A natural choice for ^^^Yb is /=^Pi, 
with a line width 7j/27r ~ 28 MHz much larger than any 
Zeeman splitting in ^ or e (thus prohibiting independent 
detection of atoms depending on their spin i / t)- Cru- 
cially, atoms in the e manifold are not in resonance with 
the imaging light and are therefore unaffected. The e^ 
atoms are subsequently brought down to the g manifold 
using, e.g.^ adiabatic passage techniques, leaving the e^ 
state unaffected. A further imaging pulse allows to detect 
those atoms, initially excited by the probe pulse. One 
might worry that a fraction of atoms from the e^ state 
could end up being transferred too, thus contaminating 
the final edge signal. Fortunately, the off-resonant ex- 
citation rate to "wrong" states will be smaller than the 
resonant rate by a factor scaling as ~ (Q/Az)'^^ with 
^z ^ ISe — 9g\B/2 a typical Zeeman splitting. Tak- 
ing for example the parameters given in [12], one has 
J/h ^ 100 Hz, and Q - 0.05J/h - 27r x 20 Hz, making 
the final contamination of ^^ by e^ negligible (~ 10~^). 



